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ABSTRACT 

When a normal shock wave moves into a still gas and strikes a blunt 
obstacle, the wave is diffracted, and a shock layer is formed on the upstream 
side of the obstacle. This shock layer adjusts with time until steady flow is 
established. 

This paper presents numerical results for the normal-shock-wave— sphere 
interaction problem for perfect- and real -gas flow fields. 

A finite difference method is used to solve the governing flow equations. 
The differencing technique which is employed is a modification to the technique 
which P. D. Lax developed for treating inviscid, unsteady one -dimensional -flow 
fields containing embedded shock waves. 

Results are presented for the unsteady behavior of the shock detachment 
distance and the stagnation-point pressure. The cases which are treated cover 
wide ranges of the ratio of specific heats and the incident shock Mach number. 
The results show that the growth of the shock- layer thickness is much slower 
than the adjustment of the stagnation-point pressure to its steady value. 
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A STUDY OF THE TRANSIENT FLOW FIELD AHEAD OF A SPHERE 
WHICH HAS BEEN STRUCK BY A NORMAL SHOCK WAVE 
By Richard W. Barnwell 
NASA Langley Research Center 

INTRODUCTION 

The transient flow which occurs when a normal shock wave moves into a 
still gas and strikes a sphere is to he investigated. This type of flow occurs 
in a shock tube which is being used as a supersonic wind tunnel during the 
period when steady flow about the model is being established. The flow field 
is assumed to be inviscid and adiabatic. Numerical results are presented for 
the time histories of the stagnation-point pressure and the shock-layer thick- 
ness along the stagnation streamline for a number of perfect-gas cases and two 
equilibrium-air cases. 

The shock-wave — blunt -body interaction problem can be approximated by 
assuming that the body is placed in a still gas and accelerated instantaneously 
to supersonic speed. Cabannes (l ) , Bausset ( 2 ), and Miles, Mirels, and Wang (3) 
have used this assumption to obtain approximate analytical expressions for the 
time history of shock-layer thickness. Experimental results have been obtained 
for the transient behavior in the flow fields of various blunt bodies which 
have been washed by normal shock waves by Davies (t ) , Bryson and Gross ( 5 ), and 
Offenhartz and Weisblatt (6) . Numerical results for diffraction of normal 
shock waves by a flat-nosed cylinder and a rectangular cylinder with its axis 
perpendicular to the flow have been obtained by Butler (j), who employs a 
numerical technique similar to the one which is used here. 
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The Numerical Method 


In the present investigation, numerical solutions are obtained to the 
initial value problem for the partial differential equations which govern 
unsteady inviscid gas flow and the boundary conditions for flow over a sphere. 
These solutions are calculated with a finite difference method which is based 
on a technique of first-order accuracy that P. D. Lax (8) developed for 
treating one-dimensional unsteady inviscid flow fields containing embedded 
shock waves. The virtue of the Lax technique is that the shock waves are not 
represented by apparent discontinuities in the flow property profiles, but 
rather by continuous profiles with steep gradients in the vicinity of the 
shocks . 

The natural coordinate system for the present problem is a spherical polar 
system with its origin located at the center of the sphere and its axis directed 
along the stagnation streamline. The governing partial differential equations 
are used in conservation form. The equations which are employed at points away 
from the axis in the present work are: 

Continuity: 

^(pr) + A.( pur ) + iL(pv) + pu + pv cot 0=0 (1) 


Radial momentum: 
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Tangential momentum: 
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Energy: 


j^-(pEr) + jp-(puHr) + ^-(pvH) + puH + pvH cot 9=0 
at or o9 


( 4 ) 
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where 


E = e + + v 2 ^ 

H = h + i^u 2 + v 2 ) 


In these equations, r and 9 are the standard polar coordinates, u and v 
are the velocity components in the r and 9 directions, respectively, p, 
p, e, and h are the pressure, density, internal energy, and enthalpy of the 
gas, and t is the time. 

Indeterminacies occur in equations (l) to (4) on the axis (9 = 0), where 
v vanishes and cot 9 is infinite. These indeterminacies can he removed with 
the aid of 1' Hospital's rule. For example, the continuity equation for 9=0 
is given hy 


(pr) + J^(pur) + 2 A( P v) + pu = 0 


The general form of the governing equations is 


SA^ SB^ SC^ 
St + Sr + S9 


+ D ± = 0 


( 5 ) 

(6) 


The Lax technique can he extended to the problem which is being treated hy 
replacing the partial derivative SAj_/St in equation (6) hy the forward finite 
difference expression 


at(( A i )^ - + ( A i )Lr,0 + 


and substituting the central difference expressions 
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and 


1 ( c \t 
2A0 \ 
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r, 0-A8 


for the partial derivatives cJBi/dr and ^C±/hB, respectively. The quan- 
tity K in the forward finite difference expression has the value sec 0 for 
the momentum equations (eqs. (2) and ( 3 ) above). For the continuity and energy 
equations (eqs. (l) and (4), respectively) K has a value of unity. This 
modification to the finite difference expressions for the time derivatives of 
the momentum equation is made in order that they will vanish for uniform flow 
parallel to the axis. The average value 
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is used to represent the undifferentiated term in equation (6) . 

When these difference expressions are substituted into the partial dif- 
ferential equations and expanded about the point (r,0,t), the following equa- 
tion results: 


dA ± , SBi 
St Sr 
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This equation differs from the general form of the partial differential equa- 
tions given by equation (6) because of the presence of the terms on the right- 
hand side. These terms cause the inviscid flow properties to be continuous 
across shock waves so that computations can be made in the vicinity of shocks 
in the same manner as elsewhere. 

The procedure which is used in this investigation consists of using the 
forward time difference given above and another time difference given by 
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on alternate time intervals. When the expansion for the first form of the time 
differences in equation ( 7 ) is replaced by that for the second form, the 
following equation results: 




Te" 
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(8) 


It has been found that the solutions which this alternating procedure yields 
are more accurate than those which are obtained when the first form of the time 
differences is used exclusively. 

It has been shown that the Courant -Friedrichs -Levy stability conditions 
are sufficient for both of the computational procedures mentioned above if 
allowance is made for the mixing of the flow properties in the vicinity of the 
smeared shocks. 


RESULTS 

Detailed Time-History Calculations 

The diffraction of the incident shock wave by the sphere and the subse- 
quent motion of the reflected shock wave for 7 = 1.4 and M s = 5*55 

= 1.70) are shown in figures 1 and 2, where the Mach numbers M s and 2^, 
are those of the incident shock wave and the uniform flow behind the incident 
shock, respectively. The nondimens ional time which is used is r = ct/ 8 ^. 


where t is the time after shock impingement, c is the initial speed of the 
reflected shock, and 800 is the steady shock detachment distance. The shock 
waves shown in the figures are positioned at the midpoints of the rapid rises 
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in the pressure profiles across the shocks. The confutation is initiated at 
r = -0.21 before shock impingement has occurred. The shock wave moves down- 
stream, strikes the sphere, and is reflected. Initially, the reflection is 
regular. When the angle of inclination of the incident shock to the surface 
has increased, Mach reflection occurs. Examples of regular and Mach reflection 
are shown in figure 1 for t = 0.088 and t = 0.37^ respectively. 

The movement of the reflected shock wave as steady flow is approached is 
shown in figure 2. The numerical results for the shock wave location at 
r = 5*28, and the experimental results of Ladenhurg, Winekler, and Van 
Voorhis ( 9 ) are in agreement near the axis where the numerical solution has 
reached the steady state. The numerical solution is still slightly transient at 
points away from the axis and in the vicinity of the shock wave. 

In figure the transient stagnation streamline pressure distributions are 
shown. The pressure drops rapidly after the initial overpressure so that the 
pressure profile within the layer approaches its steady configuration before the 
shock detachment attains its steady value. The pressure profile at r = 0.12 
shows the maximum value of the stagnation-point pressure that is encountered 
during the computation. It is seen in the figure that the maximum overpressure 
compares favorably with the pressure behind the reflected shock in the one- 
dimensional shock reflection problem. The time r = 0 corresponds to the time 
at which the incident shock should reach the body according to an exact one- 
dimensional shock computation. The numerically determined overpressure occurs 
a little later because the incident and reflected shock waves are smeared over 
several mesh spaces, and it takes several time steps for the reflection to 
occur. 
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The pressure distributions at r = 5«28 as determined by employing the 
two forms of the time differences alternately and by using the first form 
exclusively are compared in figure 3« The alternating procedure is seen to 
yield the best results. The stagnation -point thermodynamic properties are 
less in error, and the profiles are smeared less at the shock wave. It should 
be noted that the mesh spacings which are used for the two computations are 
identical; only the methods of computation differ. 

Shock Detachment Distance Histories 
Numerical results for the transient behavior of the shock detachment 
distance for several perfect-gas cases are presented in figure 4. The nondi- 
mensional quantity 8/8^ is plotted against t where the steady shock 
detachment distances are determined from the results of Van Dyke and 

Gordon (10) and Lomax and Inouye (ll). The cases for 7 = 1.1 and 7 = 1.4 
span a wide range of incident shock Mach number M s . As approaches unity 

the shock layer becomes large and difficult to compute when the steady state 
is approached. For this reason, the case 7 = 1.4 and M s = 2.51 (M® = 1.20) 
is terminated before steady flow is established. Since Mo is small for all 
values of M s when 7 is large, only one case is presented for 7 == 5 / 

Results for two computations are shown in figure 4 for the case 7 * 1.1 
and M s = 00 = 4.26). Finer mesh spacings are used for the computations 

indicated by the circles. The results show that the shock detachment histories 
are not very sensitive to the size of the mesh spacings. 

All the curves in figure 11 lie in a fairly narrow band. This indicates 
that the nondimensional time t is an effective scaling parameter. In general, 
the curves are arranged in ascending order within the band according to the 
value of Moa • 
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Several cases are treated for 7 = 1.4 to determine the asymptotic 
behavior of the shock detachment histories as M s is increased. It is seen 
in figure 4 that the histories approach the limiting case for fairly low values 
of M s ; the curves for the cases M s = JO and M s - 3-1 coincide for all 
means and purposes. In addition, the history of the shock detachment distance 
for the intermediate case 7 = 1.-4- and M s = 5*55> 'the case which is illus- 
trated in figures 1 , 2, and 3, is represented by this curve. 

Since the curves in figure 4 lie in a narrow band although 7 varies 
over a wide range, it would appear that real -gas effects should not affect the 
shock detachment distance histories appreciably. Two cases for equilibrium 
air were treated to check this observation. These cases are illustrated in 
figure 5* It is seen that the curves for these cases lie within the band of 
perfect-gas solutions shorn in figure k. 

A perfect-gas case is also shown in figure 5 . The governing parameters 7 
and M s for this case were chosen such that the Mach number and the ratio 

of enthalpy to internal energy in the region of uniform flow behind the inci- 
dent shock wave, would match those of the equilibrium case where U s = 10,000 fps, 
p Q = 10~2 atmosphere, and T 0 = 32° F. It is seen that the curve s for the 
shock detachment distance histories for these cases coincide. 

The equilibrium cases shown in figure 5 were computed with the aid of a 
curve fit of the thermodynamic properties. It can be seen from the governing 
differential equations, equations (l) to (4), that the finite difference equa- 
tions determine the density p and the internal energy e at each new grid 
point. The curve fit is used to determine the pressure p and the enthalpy h 
in terms of p and e. 
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Stagnation-Point Pressure Histories 

The numerical results for the stagnation-point pressure histories for 
several representative cases are presented in figure 6. The nondimensional 
pressure Pstag/Poo I s plotted against t. Initially, an overpressure occurs 
at the stagnation point. Then the pressure drops rapidly and approaches its 
steady value monotonically. The results of Butler ( 7 ) show that the transient 
stagnation-point pressure for a flat-nosed cylinder undershoots the steady value 
when 7 is low and M s is high. This phenomenon is not observed in the pres- 
ent results for the sphere. The maximum overpressures are in fair agreement 
with the pressures behind the reflected shocks for the one -dimensional shock 
reflection problems with similar values of M s and 7 . The numerical results 
for the steady stagnation-point pressures are in good agreement with the exact 
values. 

The stagnati on -point pressure histories for the perfect-gas case M s =9*7 
and 7 = 1.27 and the equilibrium-air case U s = 10,000 fps, p Q = 10 -2 atm, 
and T 0 = 32° F are compared in the figure. It is seen that these histories 
are similar. 

A comparison of the results in figures 7 and 5 with those in figure 6 
shows that the stagnation-point pressure approaches its steady value more 
rapidly than the shock detachment does. The pressures are within about 10 per- 
cent of the final values for r = 1. At this time, the shock detachment dis- 
tances have reached only about 50 percent of their final values. The shock 
detachments approach their final values at r = 5 or 6. The asymptotic behavior 
of the other stagnation-point thermodynamic properties is similar to that of the 
pressure. 


9 



The mesh spaeings A9 and Ar are chosen such that 1° =■ A0 < 2° and 
Ar % S ro /10 . It has been found that the shape of the curves for the stagnation 
point pressure histories is influenced by the mesh spaeings . The waviness of 
these curves is reduced if the size of the mesh spaeings is reduced. 

Approximate Analytic Solution 

A closed-form solution can be obtained for the shock detachment distance 
when Moo is large and 7 is near unity. For these conditions, it can be 
assumed that the shock is concentric and the tangential velocity gradient 
varies linearly with distance behind the shock in the stagnation region. The 
density varies only slightly with position and time. The continuity equation 
is integrated across the shock layer using the one-strip method of integral 
relations. The time -dependent detachment distance is given by 

5 - B„ (l - e" Kr ) [l - (1 - K)e- KT ] 

where 

K-l - Pb 

The quantities subscripted with b and s are evaluated at the body and the 
shock, respectively, for steady flow. 

The approximate solution for 7 = 1.1 and M s = 00 is shown in figure 4. 

CONCLUSIONS 

It is shown that the transient flow about a sphere which has been washed 
by a normal shock wave is not coupled. The stagnation-point pressure adjusts 
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to the steady value more rapidly than the shock detachment distance. In addi- 


tion, it is 
perfect -gas 


shown that the shock detachment distance histories of a number of 
cases can be scaled effectively. 
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Figure 1.- Normal shock wave- sphere interaction for 
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Figure 2.- Subsequent transient motion. 


NUMERICAL SOLUTION 
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Figure b.- Perfect gas shock detachment distance histories. 



Figure 5 *- Equilibrium air shock detachment distance histories. 
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Figure 6.- Stagnation point pressure histories. 




